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We propose a theory which describes the mechanical behaviour of magneto-sensitive elastomers (MSEs) under 
a uniform external magnetic field. We focus on the MSEs with isotropic spatial distribution of magnetic particles. 
A mechanical model is used in which magnetic particles are arranged on the sites of three regular lattices: sim- 
ple cubic, body-centered cubic and hexagonal close-packed lattices. By this we extend our previous approach 
[Ivaneyko D. et al., Macromolecular Theory and Simulations, 2011, 20, 41 1] which used only a simple cubic lat- 
tice for describing the spatial distribution of the particles. The magneto-induced deformation and the Young's 
modulus of MSEs are calculated as functions of the strength of the external magnetic field. We show that the 
magneto-mechanical behaviour of MSEs is very sensitive to the spatial distribution of the magnetic particles. 
MSEs can demonstrate either uniaxial expansion or contraction along the magnetic field and the Young's mod- 
ulus can be an increasing or decreasing function of the strength of the magnetic field depending on the spatial 
distribution of the magnetic particles. 
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1. Introduction 

Magneto-sensitive elastomers (MSEs), also known as magnetorheological elastomers, are high-tech 
materials that can change their shape and mechanical behaviour under external magnetic fields 
Nowadays, MSEs have found a wide range of industrial applications in controllable membranes, rapid- 
response interfaces designed to optimize mechanical systems and in automobile applications such as 
adaptive tuned vibration absorbers, stiffness tunable mounts and automobile suspensions ifl. 

Usually, MSEs consist of micron-sized magnetic particles dispersed within an elastomeric matrix. The 
spatial distribution of the particles inside an elastomer can be either isotropic or anisotropic (chain-like, 
plane-like), depending on the method of preparation 1%. MSEs with isotropic distribution of magnetic 
particles are synthesized by cross-linking of a polymer melt with weU-dispersed magnetic particles with- 
out any external field. To obtain an MSE with chain-like distributions of particles, one should apply a 
strong uniaxial external magnetic field to a polymer melt before and during its cross-linking [1]. Using 
a complex magnetic field with a rotating vector of magnetic strength or a strong shear flow before the 
cross-linking procedure, one can produce an MSE with plane-like distributions of particles fioll . 

The magnetostriction effect (i.e., magneto-induced deformation) and the change of mechanical mod- 
uli under external magnetic field are the most significant properties of the MSEs jllll . The deformation 
of the MSEs can be either positive (elongation) or negative (contraction) with respect to the direction 
of the applied external magnetic field. Magnetostriction and the mechanical moduli of the MSEs have 
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been investigated by experimental Illl - Eill . theoretical llTiJioll and simulation Ell-Eil studies. In theo- 
retical studies of the mechanical behaviour of the MSEs, different analytical approaches were proposed, 
which can be divided into two groups: continuum-mechanics approach and microscopic approach. In 
the continuum-mechanics approach, electromagnetic equations are coupled with appropriate mechani- 
cal deformation equations. Thus, macroscopic homogeneity of magnetic media is assumed. This approach 
predicts positive magnetostriction for such complex magnetic object as MSEs. However, the continuum- 
mechanics approach is not capable of describing a local discrete spatial distribution of particles. 

In the case of microscopic approach, magnetic particles are considered to be separated by a non- 
magnetic matrix. Dipole-dipole interaction between the particles leads to pair-wise attraction and to re- 
pulsion of the magnetic particles depending on their mutual positions. Since dipole-dipole interactions 
are very sensitive to the particle positions, spatial distribution of the particles inside the matrix strongly 
affects the sign of magnetostriction, as it was shown in the simulation (2J] and in the experiment . For 
instance, MSEs with isotropic spatial distribution of particles demonstrate an expansion along the mag- 
netic field, while MSEs with chain-like distribution of particles demonstrate a uniaxial contraction 

Recently, we have studied magnetostriction of the MSEs with isotropic and anisotropic spatial distri- 
butions of magnetic particles within the framework of a microscopic approach fid]. To describe spatial 
distribution of particles we have used a regular rectangular lattice model, which permits to consider 
"isotropic", chain-like and plane-like structures of particles. Such a regular rectangular lattice model 
predicts a negative magnetostriction of MSE for all distributions of particles. Our prediction for nega- 
tive magnetostriction is in agreement with experimental works (isl |27|, |28|] for chain-like distribution. 
However, predictions of the rectangular lattice, which degenerates into cubic lattice for isotropic spatial 
distribution of particles, contradict the experimental data for magneto-induced deformation of the MSEs 
with an isotropic distribution of magnetic particles L20,,28J. The regular rectangular lattice is probably 
not capable of sufficiently describing the isotropic distribution. The reason is that any regular lattice is 
intrinsically anisotropic, since one always finds the lattice directions with different regular distances be- 
tween the lattice sites. However, in spite of this feature, some regular lattices can possibly be used to 
reasonably well model the mechanical behaviour of the isotropic MSE. Therefore, in the present study 
we consider different lattices to describe isotropic distributions of magnetic particles in an MSE: simple 
cubic, body-centered cubic and hexagonal close-packed lattices. For these three different lattice models 
we examine magnetostriction and Young's modulus of the MSE in the presence of an external magnetic 
field. We construct the free energy which consists of elastic and magnetic energies. The Neo-Hooke law is 
used to describe entropic non-linear elasticity of polymer chains. 



2. Free energy of MSE 

2.1. General equations 

In this section we derive an expression for the free energy of an MSE. Analysis of the free energy 
will permit us to study the mechanical properties of the MSE. Typical MSEs are characterized by a strong 
elastic coupling between the particles and the matrix. This coupling is due to the fact that magnetic parti- 
cles are much larger than the mesh size of the polymer network. Since our paper is devoted to the static 
mechanical properties of MSEs, we consider only average displacements of particles coupled to the ma- 
trix. Under these assumptions the free energy of the deformed MSE under external magnetic field can be 
written in the form: 

F=Fei + Fm. (2.1) 

The first part Fei is the elastic energy of a deformed MSE due to the entropic elasticity of polymer 
chains. In the present work we extend a formalism presented in reference |26] for linear deformations 
of MSEs to non-linear deformations. To calculate the elastic part of the free energy _Fei of an MSE under 
finite deformation as a function of the strain we use the approach of a continuous medium. This approach 
means that a sample is divided into representative volumes of an elastic matrix which contain a large 
number of particles. Under this assumption, the free energy Fei of an MSE under non-linear deformations 
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can be expressed through the Neo-Hooke law that can be written in a general form 12911 : 

^^ei=Y(/i-3), (2.2) 

where the material parameter Go is the effective shear modulus of an MSE, Ii is the first scalar invariant 
of the Finger strain tensor B: h = trB. The value of Go includes contributions of different possible effects 
into the elastic energy appearing under elongation of a sample: reinforcement of an elastic matrix by 
hard particles due to non-uniform deformation of the matrix between the particles (strain amplification), 
possible adhesion of a polymer matrix onto the surfaces of hard particles, etc. However, we do not discuss 
here how the value of Go depends on these effects, since this task is a special problem in the theory of 
elasticity for isotropic reinforced rubbers L301 . Instead we use Go as a phenomenological parameter of the 
theory assuming that it can be extracted from experimental data for elasticity of an MSE in the absence of 
a magnetic field. In our theory we assume an incompressible elastomeric matrix with Poisson's coefficient 
V = 1/2. Due to the axial symmetry with respect to the external magnetic field H, the MSE will provide an 
uniaxial elongation along H. In this case, the first scalar invariant 7i has the form: 

h^Z^l (2.3) 

i 

where Ajc, ^y, are the elongation ratios for the deformation of an elastomer in the three principal 
directions. 

The mechanical response of an elastomer to the magnetic field is characterized by the value of the 
strain e = AZ/Z, where AZ and Z are the elongation and the original size, respectively, of an elastomer 
along the direction of the magnetic field (x-axis). The condition of constant volume for elastomers (illliill 
permits us to relate the elongation ratios Xx, Ay, Az as follows: 



A;c = l + £, Ay = Az = 1/ Vl + £. (2.4) 

Within the framework of Neo-Hooke approximation, the elastic energy for elastomer as a function of 
£ can be written as: 



Pel^ — 
b 



(l + £)2 + - 3 

l + £ 



2 



r-2 



-^ + 0{E^), (2.5) 



where £0 is the Young's modulus of a filled matrix: £0 = 3Go (si]. In the present theory, the value of the 
Young's modulus £0 of the filled elastomer will be considered as a parameter in calculations. 

The second part of the free energy in equation F^n arises from the potential energy of magnetic 
particles placed in an external magnetic field. Application of a magnetic field induces an average mag- 
netic moment in each particle. The values of the induced magnetic moments in the magnetic particles 
depend on the material of the particles. Magnetic particles that are usually used in preparing MSEs have 
micron-sizes and a multi-domain magnetic structure. Nevertheless, they are found to show very narrow 
hysteresis cycles which indicates a soft magnetic behaviour. In this case, the dependence of magnitude 
of the particle magnetization M on the field strength H can be described in a good approximation by the 
FrohUch-Kennely equation jillisll 

M= , (2.6) 

Ms + (Aiini-l)|H| 

where Ms is the saturation magnetization and /iini is the magnetic permeability of the particles. The 
values of the saturation magnetization and magnetic permeability for carbonyl iron particles are well 
established in experiment and are equal to Ms =; 1582 kA/m and fiini ~ 21.5 with the average diameter of 
particles of 470+ 180 nm (H]. Similar values were obtained for particles of the size of 2iim: Ms=1990kA/m 
and ;L<ini - 132 fisi]. Equation <2.6> can be rewritten in terms of dimensionless parameters: reduced mag- 
netic field h - (/iini - 1)H/Ms and reduced magnetization MIM^ in the following form: 

M h 

— = . (2.7) 

Ms \ + \h\ 

The dependence of M/M^ on h is given in figure [l] Equation I l2.7t is used below for calculations of me- 
chanical characteristics of MSEs in an external magnetic field. 
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(H^„-l)///M, 

Figure 1. Reduced magnetization Ml Ms as a function of a reduced magnetic field h = ifimi - l]HI Ms- 
Magnitude of a reduced magnetization increases with an increasing magnetic field and tends to satura- 
tion magnetization Ms, when H ^ oo. 



The interaction energy of magnetic particles in an external magnetic field consists of two contri- 
butions: the dipole-dipole interaction energy between the particles and the dipole-field interaction en- 
ergy (37, 38]. However, under our assumptions given by equation I l2.6t . the dipole-field interaction energy 
- (M • H) is independent of the strain e, since the bulk magnetization of the sample M is assumed to be 
independent of strain £. Therefore, the dipole-field interaction energy does not provide a contribution 
to the mechanical characteristics and it can be excluded from our considerations. Thus, in our approach 
the magnetic energy Fm includes only dipole-dipole interaction energy per unit volume, which can be 
written as: 

' 3 (rh; ■ R,-^- ) {ihj ■ R,-; ) (m,- ■ m ^ ) 



V 4n ^ 



(2.8) 



where |Uo is permeability of the vacuum, V is the volume of the sample, /Xr is the relative permeability 
of the medium. In the present work we consider an elastomeric matrix to be non-magnetic. Therefore, 
ever5rwhere below we set /ir - 1. Here m,- and rhj are dipole moments of i-th and 7-th magnetic particles, 
R,y is the radius vector that joins the i-th and j-th magnetic particles. The radius vectors Rj-y depend 
on the macroscopic deformation £. To relate Rtj with £ we use an approximation of affinity of deforma- 
tion (illliil], as well as incompressibility of elastomers: 

(2.9) 
(2.10) 
(2.11) 

where {Ri])( and are the components of vectors that join two magnetic particles after and before 
deformation, respectively, - x,y,z). Summation in equation <2.8) is performed over all pairs of parti- 
cles. To simplify the calculations, we consider regular spatial distributions of particles. In particular, we 
use three different lattices in order to study the effects of particle distribution on the magneto-mechanical 
properties of the MSEs in a homogeneous magnetic field. 



ij)x - 
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2.2. "Isotropic" distributions of magnetic particles: different lattice models 

In contrast to our previous study (2^ . which used only a cubic lattice to describe the isotropic distri- 
bution of magnetic particles inside the MSE, we consider here different lattice models (see figure [2). In 
these models, it is assumed that the magnetic particles in an isotropic non-deformed MSE are located at 
the sites of simple cubic (SC), body-centered cubic (BCC) and hexagonal close-packed (HCP) lattices. Let a 
be the edge length in the three lattices [figure |2](b)]. 
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Figure 2. Schematic illustration of an MSE with isotropic spatial distribution of particles (a) and the lattice 
models (simple cubic, body-centered cubic and hexagonal lattices) for modelling the mechanical proper- 
ties of an MSE (b). 



For simplicity we assume that all particles are the same and have a spherical form; r is the radius of 
particles. The value of r characterizes the average size of particles in a real elastomer. Then, the volume 
fraction, (p, of the particles distributed in the matrix is given by: 

= r-?. (2-12) 
where vq - \nr^ is the volume of a particle. The factor 7 depends on the type of the lattice: 

rsc = l, rBCC = 2, 7hcp = 2V2. (2.13) 

Due to the symmetries of the infinite lattices to be used (translational symmetry and the axial symme- 
try around the vector H), the induced magnetic dipoles m, and my are directed along the external field H 
[x-axis in figure[2](b)] and their absolute values are identical m,- - mj - vqM, with M being the magneti- 
zation of each particle. Besides, a regular arrangement of the magnetic particles on the sites of the lattices 
makes it possible to simplify the summation in the equation <2.8) . To perform summation over indexes / 
and j we choose j-th particle and make summation over index /. Due to the translational symmetry for 
infinite lattices, the contribution to the magnetic energy of the dipole-dipole interaction of a given 7-th 
particle with other particles around it does not depend on the number j. Thus, the double sum over pairs 
of indexes ; and j is reduced to a simple sum over index / multiplied by N, where N is the number of all 
particles. Then, equation l l2.8t can be rewritten in the form: 

where c - N/V is the concentration of particles. Here, for convenience we calculate the radius vectors 
from the origin of a lattice and simplify the notation as R,y = R,. We introduce the parameter uq: 

Mo=^=^, (2.15) 
An 

that defines the characteristic density of energy of magnetic interaction. For Ms ~ 2 x 10^ A/m we obtain 
Mo = 4 X 10^ Pa. Below, we will show that mechanical behaviour of an MSE in the magnetic field is deter- 
mined by a dimensionless parameter EqIuq, i.e., by the ratio between characteristic values of the elastic 
and magnetic energies. 



3(R,)^-|R,f 



(2.14) 
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The radius vector R; in equation 1 12.14) is related with the radius vector R; in a non-deformed MSE 
by equations <2.9) - t2A]l . For simple cubic, body-centered cubic and hexagonal close-packed lattices, the 
value Rj can be presented in the form: 

R-=a-r;. (2.16) 

The dimensionless vector runs now over all sites of the infinite lattice with the edge length a - 1 
except f J - 0. Using the Frbhlich-Kennely equation for reduced magnetization Ml Ms, we express 
the magnetic energy Fm as a function of a reduced magnetic field h and the strain e: 



h 



l + \h\ 



/(e), 



where the dimensionless function /(e) has the following form: 

1 3 „ 2(l + e)3(r,-)2-(ri)2-(r,02 

/(£) = --(l + e)5 X ^ 



(2.17) 



(2.18) 



We recall that the vector f runs over the sites of different lattices with a - 1. The value of function f{£) 
depends on the type of the lattice, since the components of the vectors f ,• depend on the lattice. The vector 
for three lattices can be expressed in the following form: 



?/ = d-ia + b-ib + c-ic 



(2.19) 



where ia, it and ic are the integer numbers, a, b and c are the basis vectors. The basis vectors are 
presented in figure[3]and have their coordinates: 






X 






















/ b 






Figure 3. (Color online) Basis vectors a, b and c introduced for the simple cubic, body-centered cubic and 
hexagonal close-packed lattices with the unit edge length, |b| = 1. 



Simple cubic lattice 

a= (1,0,0), S= (0,1,0), c= (0,0,1) 
Body-centered cubic lattice 

7 ^ fl 1 1 

fl= (1,0,0), fo= (0,1,0), c= 



(2.20) 



(2.21) 



• Hexagonal close-packed lattice 



2 1 1 



, V 3 2 



fo= (0,1,0), 



( 1 v^^ 
0,-,— 
2 2 



(2.22) 



For different lattices, the function /(c) was analyzed numerically. Below we consider the magne- 
tostriction and the Young's modulus of the MSE in the presence of an external magnetic field depending 
on the type of the lattice model. 
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2.3. Free energy as a function of deformation 



Using equation for elastic energy and equation <2.17> for magnetic energy, we obtain free energy 
of the MSE in the following form: 



Eq 
6 



(l + £)^ + 



l + £ 



h 



l + \h\ 



/(£), 



(2.23) 



where the values of the function /(e) vary for the simple cubic, body-centered cubic and hexagonal 
close-packed lattices. We calculate the reduced free energy, FIuq, of an isotropic MSE as a function 
of the strain e for the three types of initial lattices at fixed values (f) - 0.05 and EqIuq - 2.5 and at 
different values of the reduced magnetic field h - 0,0.25,0.70,1.5,4,40 that correspond to the values 
M/Ms - 0,0.2,0.4,0.6,0.8, 1.0, respectively. One can see in figure|4]that application of the magnetic field 
leads to the shift of the minimum of free energy. The minimum of free energy determines the equilib- 
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Figure 4. (Color online) Reduced free energy, Fl mq. of an isotropic MSE as a function of the strain e calcu- 
lated for three types of initial lattices: simple cubic lattice (a), body-centered cubic lattice (b) and hexago- 
nal lattice (c) at different values of the reduced magnetic field h = (/ijni - l)HIMs. 



rium elongation of an MSE, Cgq- One can see that the value and the sign of the Ceq depend on the lattice: 
£eq < for the simple cubic and hexagonal close-packed lattices and eeq > for the body-centered cubic 
lattice. This means that both simple cubic and hexagonal close-packed lattices predict a contraction of an 
isotropic MSE along the magnetic field H (eeq < 0), whereas the body-centered cubic lattice predicts an 
expansion of an isotropic MSE along the magnetic field H (Ceq > 0). 

The analysis of Ceq as a function of the reduced magnetic field h, the volume fraction cp and parameter 
Eq/uq for different lattices are given below. 



3. Mechanical behaviour of MSEs in a homogeneous magnetic field 
3.1. IVIagnetostriction of the IVISE 

In this section we study the magnetostriction effect in the MSE with isotropic spatial distribution 
of magnetic particles. We calculate the equilibrium elongation of the MSE from the minimum of free 
energy _F: 



dF 

ae 



= 0. 



(3.1) 
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Figure 5. (Color online) Dependence of the equilibrium elongation Ceq on the reduced magnetic field 
h = (^ijjj - l]HIMs at different volume fractions (p, calculated for three types of initial lattices: simple 
cubic lattice (a), body-centered cubic lattice (b) and hexagonal close-packed lattice (c). 



Using equations <2.23> and <2.18> . the equation 13. 1> for the equilibrium elongation Eeq can be rewritten 
as follows: 



Eo 
3 



1 + e, 



1 



eq- 



(l + eeq)2 



h 



l + \h\j 7 



eq 



^ 12(1 + EeqfirOt - mri)la + Eeqf [infy + inf.] + 3 [ [nfy + (r;)^] ^ ^ 

2[(l + eeq)3(r;)2 + (ri)2 + (r;)l]^ 



(3.2) 



Dividing both the left- and right-hand sides of equation i3.2\ by the factor uq, one can see that the equi- 
librium elongation Ceq depends on the elastic modulus £0 and on the magnetic parameter uq through 
their dimensionless ratio EqIuq. The last equation we solve numerically at varied values of the reduced 
magnetic field h. Figure [5]shows dependences of the equilibrium elongation £eq on the reduced magnetic 
field h at the values of parameter EqIuq - 2.5 and volume fraction (p - 0,0.01,0.05 and 0.1, calculated 
for the three types of lattice models. Figure [6] shows dependences of the equilibrium elongation feq on 
the reduced magnetic field h at the values of volume fraction (f) - 0.05 and parameters EqIuq- 1.0,2.5,5 
and 10, calculated for the three types of lattice models. One can see that for simple cubic and hexagonal 
close-packed lattices, an MSE is uniaxially contracted along the direction of the external magnetic field, 
Eeq < [figure[5](a), (c) and[6](a), (c)], while for the body-centered cubic lattice, an MSE uniaxially expands 
along the direction of the external magnetic field, Eeq > [figure[5](b) and[6](b)]. 

Different signs of magnetostriction of MSEs with different spatial "isotropic" distributions of particles 
(simple cubic, body-centered cubic and hexagonal close-packed lattices) can be explained by the fact that 
the mechanical behaviour of MSEs is determined mainly by the mutual attraction and repulsion of the 
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Figure 6. (Color online) Same as figure [s] but at different values of the parameter EqIuq and at fixed 
volume fraction ip = 0.05. 
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nearest particles as it is illustrated in figure [T] In a simple cubic lattice, the nearest particles to a given 
one (B) are the particles A and C as presented in figure[7] Such configurations lead to the contraction of an 
MSE along the magnetic field which is in accord with our results presented in figures[5](a) and[6](a). In the 
body-centered cubic and hexagonal close-packed lattices there exist some extra particles whose positions 
are determined by the angle 6 (see figure[7). Depending on the value of 6, either contraction or elongation 
of an MSE is possible. One can show that interaction between the particles B and D results in the expan- 
sion at 32° < < 72° and in contraction at < < 32°, 72° <B< 90°. Therefore, the body-centered cubic 
lattice predicts an expansion since 0bcc = 54.7° . For a hexagonal close-packed lattice, the nearest particles 

which lie at the angles 0hcp - 35° provide an expan- 
sion and the nearest particles which lie at 0hcp - 90° 
provide a contraction. We have shown that the inter- 
play between two different contributions of the near- 
est particles leads to the contraction. Thus, the sign of 
magnetostriction strongly depends on the local spatial 
distribution of magnetic particles. 

Furthermore, one can see in figure [s] that the in- 
crease of the volume fraction (p results in an increase 
of the equilibrium elongation Iceql, when h is fixed. Fig- 
ure [6] shows that the increase of the parameter Eq/ uq 
results in a decrease of matrix deformation |eeqL when 
h is fixed. These results are explained by the fact that 
the relative contribution of magnetic interaction be- 
comes larger at higher values of (p and at lower values 

of Eq/Uq. 




H 



repulsion 



Figure 7. Attraction and repulsion of magnetic 
particles in an MSE depending on their mutual 
positions. 



3.2. Young's modulus of the MSE 

In the case of a tensile deformation, we consider such geometry when the mechanical force is applied 
along the external magnetic field H, i.e., along the x-axis in our case (see figure[2). The Young's modulus 
E can be obtained as the second derivative of free energy with respect to e, where e - Ceq + Se: 



(3.3) 



which yields 
1 + 



Eq 
3 



(l + eeq)3 



l + \h\j 4v/TT^,-f 



E 



32(r,-)^(l + eeq)^- 



imrOl (ri);+(ri)^ (l + £eq)°+90(r;)^ (ri)t+(r,-)^ (l + £eq) 



(3.4) 



One can see that the ratio E/Eq depends on the parameters £0 and uq through their dimensionless ratio 

EqI Uq. 

We have numerically calculated £ as a function of the reduced magnetic field h using equation <3.4> . 
Figure[8]shows the dependence of the Young's modulus E on the reduced magnetic field h at the values of 
parameter Eq/uq - 2.5 and volume fraction cp - 0,0.01,0.05 and 0.1, calculated for three types of initial 
lattices: simple cubic, body-centered cubic and hexagonal close-packed lattices. Figure[9]is the same as 
figure|8]but at a fixed value of the volume fraction cp - 0.05 and at varied values of EqIuq- 1.0,2.5,5 and 
10. One can see that the Young's modulus E decreases for the simple cubic [figure[8](a) and[9](a)] and the 
body-centered cubic [figure |8] (b) and[9](b)] lattices and increases for the hexagonal close-packed lattice 
[figure[8](c) and[9](c)]. Different behaviour of £ as a function of a reduced magnetic field h for the three 
lattices is caused by complex combinations of interactions between particles of different mutual positions 
(see figure[7). 
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Figure 8. (Color online) Dependence of the Young's modulus E on the reduced magnetic field h = (^inj - 
I) HI Ms at different volume fractions (p, calculated for three types of initial lattices: simple cubic lattice 
(a), body-centered cubic lattice (b) and hexagonal lattice (c). 
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Figure 9. (Color online) Same as figure [8] but at different values of the parameter EqIuq and at fixed 
volume fraction ip = 0.05. 
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Furthermore, it can be seen in figure |8]that an increase of the volume fraction (p leads to an increase 
of absolute values of the change of the modulus \E- Ea\ for all distributions at fixed h. Figure [9] shows 
that an increase of the parameter EqIuq results in a decrease of absolute values of the change of the 
modulus \E- Eo\ for all distributions at fixed h. These results are explained by the fact that the relative 
contribution of the magnetic energy to the modulus increases at increasing values of the volume fraction 
(p and decreases at increasing values of the parameter EqIuq. 



4. Conclusion 

In this paper we have studied the mechanical properties of magneto-sensitive elastomers with isotro- 
pic distribution of the magnetic particles in an external magnetic field. We have used a model in which 
magnetic particles are located at the sites of regular lattices. Different types of the lattice models have 
been considered: simple cubic, body-centered cubic and hexagonal close-packed lattices. Such lattice dis- 
tributions of particles can be prepared in experiment as was shown in 1 9] . We show that magneto-induced 
deformation and the Young's modulus of the MSE strongly depend on spatial distribution of magnetic 
particles. A simple cubic lattice model predicts a contraction of the MSE with isotropic distribution of 
magnetic particles along the direction of a homogeneous magnetic field. It predicts that the Young's mod- 
ulus decreases with an increase of the magnetic field; the same result is obtained for the body-centered 
cubic model. However, in contrast to the simple cubic lattice, the body-centered cubic lattice model pro- 
vides an expansion of the MSE along the direction of a magnetic filed. The hexagonal lattice model shows 
that MSEs with isotropic particle distribution contract along the external magnetic field, while its Young's 



33601-10 



Magneto-sensitive elastomers 



modulus increases. These findings may explain different signs of magnetostriction observed in experi- 
ments with differently prepared MSEs I2sll . 
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Bniii/iB po3noAi.ny nacTi/iHOK Ha MexaHiMHi B/iacTi/iBOCTi 
MarHiTHOHyT;ii/iBi/ix eyiacTOMepiB b OAHopiAHOMy MarHixHOMy 
no;ii 

/],. iBaHeMKcl^, B.n. ToineBiKO^E] m. Ca(j)'aHHiKOB^, f. raMHpi)^^ 

iHCTWTyr MaTepia/iosHaBCTBa, flpeBfleHCbKwPi TexHi4HHw VHiBepci/iTeT, 
Byn. re/ibMro;ibL;a, 7, 01069 flpesfleH, HiMeH4HHa 

flpesfleHCbKWM iHCTMTyT no/iiMepHHX floc;iifl>KeHb iM. /lati6HiL^a, ayn. fore, 6, 01069 flpesflen, HiMeH4i/iHa 
iHCTwryr BucoKOMO/ieKy/i^pHnx cno/iyK, PociticbKa AKafleMia HayK, 
Bo/ibi±ioCi np., 31, 199004 CaHKT-neTep6ypr, Pocia 

Mn nponoHyeMO reopiio, SKa onucye MexaHinny noBefliHKy MarHiTHOHyr/inBnx e/iacTOMepiB (MHE) 3 isoipon- 

HWM pOSnOfli/IOM 4aCTWH0K B OflHOpiflHOMy MatHiTHOMy nO/li. BWKOpi/ICTOByiOTbCfl Tpn MexaHii-lHi MOfle/li, B 

51KHX MarHiiHi HacTMHKH HOMimeHO y Bys/iax BiflnoBiflHO ipbox nepioflW4HHx I'paTOK: npocToT Ky6i4H0i, o6'eMO- 
L4eHTpoBaHoT Ky6iHH0i'Ta mi/ibHO BnaKOBaHoT reKcai'OHa/ibHoT. U,v\m mm posBUBaeMO Hai±i nonepeflhiw niflxifl 
[Ivaneyko D. et al., Macromolecular Theory and Simulations, 2011, 20, 411], b HKOMy mh BHKopMCTa;in nvwue 
npocry Ky6iHHy fpajny fl/ia onncy npocTopoBoro po3nofli;iy nacTHHOK. /l,ec|)opMaL;ifl, ^Ka BMHMKae nifl bh/ih- 
BOM MarHiTHoro no/ifl, ra MOfly/ib tOnra b MHE poBpaxoBani fl/ia piBHUx HanpyxeHOcreii MarHiiHoro nons^. Mm 
noKasa/iM, mo MarHiTHOMexaHiHHa noBefliHKa MHE CM/ibHO 3a/ie>KMTb Bifl npociopoBoro po3MimeHH« MarHir- 
HMX 1-iacTMHOK. MHE fleMOHcrpye oflHOBicHe BMflOBxeHHfl a6o ctmck BSflOBX MarHiiHoro nonn, a MOfly/ib tOHra 
spocrae a6o cnaflae npM 36i;ibLueHHi HanpyxeHocii MarHiiHoro no/ia 3a;iexH0 Bifl npocropoBoro po3nofli;iy 

MarHiTHMX 4aCTMH0K. 

K/hohobI c;iOBa: MarniTHO^yTJiMBi enacTOMepn, MexaHinni BJiacTusocTi, MOffynb, MO^ejiKiBaHHa 
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